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Abstract 

The microscopic description of spin-Peierls substances is discussed. Particular attention is paid to the ordered 
(dimerised and incommensurably modulated) phases. Important points are the adiabatic and the antiadiabatic 
approach, generic soliton forms, elastic and magnetic interchain couplings. The wealth and the accuracy of exper- 
imental information collected for the first inorganic spin-Peierls substance CuGeOs motivates this work. 



Keywords: non-adiabaticity; soliton form; local renormalisation; interchain coupling 



The coupling of ID magnetic and 3D elastic 
degrees of freedom in spin-Peierls systems leads 
to magnetically driven lattice distortions at low 
T (see e.g. [1,2]). Without magnetic field the spin 
chains dimerise (D phase). For sufficiently large 
magnetic field or by defects the homogeneous 
dimerisation is modulated. Solitons occur: zeros 
of the dimerisation at which S = 1/2 spinous are 
localised [3,4]. The high-magnetic field, incom- 
mensurate (I) phase is characterised by an array 
of these solitons. 

Since we deal with coupled magnetic and elastic 
degrees of freedom the microscopic Hamiltonian is 
H — Hs + HsB + Hb- The purely magnetic part is 
given by 



treated as free bosons 



(1) 



where i is the site index within a given chain j. 
The frustrating next nearest neighbour coupling 
a J is included since there is evidence that this term 
is important in CuGcOs [5-7]. The phonons are 



(2) 



The Peierls phenomenon stems from the interac- 
tion HsB ■ Since the atomic displacements are small 
we assume that they influence the magnetic ex- 
change in linear order 



HsB=J2Mbt + b- 



with 



k 



k-q 



(3) 
(4) 



Unfortunately, this Hamiltonian is too compli- 
cated. So one argues that the ID magnetic subsys- 
tem strongly favours a distortion at wave vector tt 
which is underlined by the dimerisation below the 
transition temperature Tgp and by the diverging 
dimer-dimer correlation at tt for T — > 0. Hence 
those phonons matter most which are in the vicin- 
ity of the actual lattice distortion at for T < 
Tgp. It appears permissible to take their energy 
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a;(q'^) as the energy of all the phonons neglecting 

the phonon dispersion. This does not imply that 
the 3D character is lost since in generic lattices 
with basis even the local phonons are extended 
objects (for the polarisation vectors of CuGeOs 
see [8-10]). This fact leads to the important elastic 
coupling between adjacent chains so that the spin- 
Peierls transition can take place at T > 0. A purely 
ID model would not display a phase transition at 
finite T. By the above simplifications wo arrive at 
Hb = '^J2ijK,j^i,j- The spin-phonon coupling 
in real space reads Hsb = 9J2i,j ^ijiKj + Kj)- 
The operators Aij contain spin operators not too 
far from the point indexed by e.g. Aij = 

9{^i,j^i+l,j ~ Si,jSi_i,j + / X^-|-(Si+l,j±lSj,j±i — 

Sjj±iS,_i.j±i)) [11]. The parameter /(|/| < 1) 
measures how much a phonon on chain j infiuences 
also its adjacent chains j ± 1. 

Now two routes can be followed. The first, more 
traditional, route consists in treating the phonons 
adiabatically, i.e. the distortions are considered as 
static 



=: 6i,,J/i2g) 



(5) 



in the ordered phases. Correspondingly, there is no 
influence of Hsb above Tsp, see e.g. Ref. [12]. The 

phonon propagators arc rcnormaliscd in RPA by 
the magnetic response in the dimer-dimer channel 
[13]. The transition is signalled by the vanishing of 
the renormalised phonon frequency, the so-called 
phonon softening. This description applies if the 
phonons are the slow degrees of freedom [13] which 
are renormalised by the fast magnetic degrees of 
freedom. An obvious condition for the validity of 
this scenario is w < J. 

It was shown by numeric investigation that the 
more restrictive condition a; < A applies where 
A is the gap in the D phase [14]. The gap A is a 
measure for the effect of the spin-phonon coupling, 
i.e. A indicates up to which energy the magnetic 
and the phononic degrees of freedom are decisively 
altered by the coupling. Hence the phonon cannot 
become soft if w > A holds. This view, suggested in 
particular for the spin-Peierls substance CuGeOa 
[11], is corroborated by field theoretic results [15]. 



The second route consists in the antiadiabatic 
approach which considers the spin system as slow 
and the phononic system as fast. Then it is rea- 
sonable to eliminate the phonons in favour of an 
effective spin model. Many works took actually 
this approach [16,17] in leading order in J/uj. Cou- 
plings dependent on T arise in the next-leading 
order [11]. The transition into the ordered phase 
docs not occur due to the softening of a phonon 
but due to the tendency of the effective spin model 
towards dimerisation. Considering a single chain 
phonon-induced frustration above its critical value 
ctcff > Q!c ~ 0.241167(5) [18] drives the dimerisa- 
tion [11,19,17]. Hence a finite spin-phonon interac- 
tion is necessary to achieve the spontaneous sym- 
metry breaking [20-22] . 

In Sect. 1 we compare the adiabatic and the an- 
tiadiabatic route to the spin-Peierls transition. In 
Sect. 2 the I phase will be discussed in detail. A 
summary concludes this article. 



1. Mean-field, nonadiabaticity and phasons 

Most adiabatic approaches use the Hamiltonian 

K 



h=jj: 



(6) 



with some spring constant K (see e.g. [23,24]); the 
chain index has been dropped. The elastic inter- 
chain coupling does not appear since its main ef- 
fect can be absorbed in the ID model as long as 
the modulations on neighbouring chains occur in 
phase. An important exception to this rule are soli- 
tons induced by impurities [25,26,4]. The local dis- 
tortions Si in Eq. (6) are subject to the constraint 
Si = 0. The ground state energy (H) is min- 
imised with respect to the Si leading to 



KSi = -(Si+iSj) + (Sj+iSi) 



(7) 



where the overline stands for the sample average 
to ensure the constraint. 

For 5tot = 0, i.e. in the D phase, the solution 
to the minimisation of Eq. (6) reads Si = (— l)*(5o 



2 



where the amplitude is proportional to the ex- 
pectation value of the dimerisation operator Hd = 
Z.i-iySi+iSi, i.e. 60 = {Ht,)/K (see e.g. [27]). 
For finite magnetisation solitons occur. A generic 
solution for Sz = I is displayed in Fig. 1 [3]. 




Fig. 1. Upper panel: local distortions. Symbols: 
DMRG-result at K = 18J, a = 0.35; solid line from (18) 
with S = 0.014, fed = 0.959, = 10.5. Lower panel: lo- 
cal magnetisations; solid line from (17) with W = 0.21, 
R = 5.0, km = 0.992, = 7.9. 



What does one obtain in the antiadiabatic limit? 

For simplicity we assume Ai_j = g{: S^+i^S^j : 
: Sj+i j±iSij±i :) where the colons stand 
for the operator reduced by its expectation value (: 
X := X — {X)). This is to prevent "trivial" rcnor- 
malisation of the bare coupling constants. Apply- 
ing a unitary transformation eliminating ifsB one 
obtains in the two leading orders additional terms 
in the spin Hamiltonian AHx = 0{g'^/uj) and 
AHy = 0(<?2j/a;2) [11] 

^^^ = i^°*K^)E[4.'LA..] . (9) 



where L is the Liouville operator, i.e. the commu- 
tation with Hs- In AHy the interchain terms in 
Aij do not generate qualitatively new terms. In 
AHx, however, the interchain terms lead to terms 

Hint = — ^-^y^'^»+i'J^»'.J' ••Si-i-ij-i-iS»>j-M -(10) 

U) — 7 



coupling adjacent chains. It is obvious that via Hi^t 

the dimerisation on one chain favours dimerisation 
on the adjacent chain. So in the antiadiabatic de- 
scription the part (10) is responsible for the phase 
transition at T > since it creates the coherencx; 
of the dimerisation pattern throughout the whole 
sample. 

It is often argued that the adiabatic phonon de- 
scription works very well, see e.g. [12,10,15]. The 
pieces of evidence, Ginzburg criterion, suscepti- 
bility and entropy fits, do indeed favour a mean- 
field description. They do not, however, require 
an adiabatic phonon description. In view of cj ^ 
A in CuGeOa [11], we advocate the antiadiabatic 
phonon description. This is combined with a chain- 
mean field treatment (CMF) since Tgp/J » 0.1 is 
small. 

Within each chain an effective Hamiltonian as 
derived previously [11] is used. Between the chains 
the contribution (10) is treated in mean-field. Let 

us abbreviate 



6i := - 



U)J 



(11) 



Then the 5i can be interpreted as local distortions 
as in the adiabatic description, cf. (6). Defining 



UjJ 



45V 



(12) 



the equivalence is pushed one step further. Insert- 
ing Eq. (12) in Eq. (11) and treating both ncighour- 
ing chains as equal yields in CMF the same equa- 
tion (7) as in the adiabatic treatment. This is true 
in the D phase and in the I phase as long as the 
modulation stays close to tt. 

The "good" equivalence between the adiabatic 
treatment and the antiadiabatic treatment with 
chain-mean-field provides an explanation for the 
success of the former in spite of the high phonon 
energies. But it is not the adiabatic treatment in 
its own right which works so well. There are two 
major differences between the adiabatic treat- 
ment and the seemingly equivalent antiadiabatic 
treatment, (i) The in-chain couplings are partly 
phonon-induced and hence they depend on T 
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[11,28]. (ii) In the adiabatic treatment K is related 

to in-chain properties: K = a;J/(2g^) (obtained 
inserting Eq. (5) in Hsb and Hb). The constant 
Kgs in Eq. (12) refers to an inter-chain property 
since the spatial extent of the phonons as mea- 
sured by / enters decisively. 

What are the points in favour of the antiadia- 
batic approach for CuGeOa? Besides the relation 
between A « 25K and w « 150K [11] and the 
T-dependent couplings [28] there is the observa- 
tion that the zero-point motion of lattice exceeds 
the static distortions considerably [10]. In partic- 
ular in the I phase, this zero-point motion has an 
important effect on the alternating local magneti- 
sation [29,3]. On the basis of results as in Fig. 1 
NMR peaks can be understood [30]. In adiabatic 
calculations the amplitudes of the alternating com- 
ponent, however, is found to be much larger (fac- 
tor 4 to 6) than in experiment. Theory and ex- 
periment can be reconciled if the zero-point mo- 
tion of the "phonons" of the soliton lattice, the so- 
called phasons, are taken into account. This mo- 
tion leads to a reduction of the alternating compo- 
nent [3] which can be mimicked by an average of 
the magnetisation of adjacent sites [29]. The pha- 
sons are connected to the spontaneous breaking 
of the quasi-continuous symmetry of translating 
the modulations along the chains. Such transla- 
tions require that the distortions are also dynamic. 
So the phasonic zero-point motion provides strong 
evidence for the nonadiabatic character of the or- 
dered phases in CuGeOa. 



2. Incommensurate phase 

Turning to the solitons in the I phase we first 
state the standard theory [31]. The original spin 
Hamiltonian (6) is mapped by a Jordan- Wigner 
transformation [32] to a fermionic model which in 
turn can be reduced to an effective bosonic contin- 
uum model at low energies [31,33] 



+ J (-5(x)cos(20) + — (52(x))dx , (13) 

where S(x) is the alternating component of the 
distortion. Nakano and Fukuyama used the self- 
consistent harmonic approximation (f) ^ (pel + ^fl 
distinguishing a classical number (j)ci and a fluc- 
tuating operator part (pa with {(pn) = 0. Treat- 
ing the fluctuation in harmonic approximation 
leads to cos(20) e"^'^ cos(20ci)(l - 2(^^) with 
cr {4'fi{x)) motivated by the WKB approach 
by Dashen et al. [34]. Varying the ground state 
expectation value of (13) with respect to d{x) and 
to (pel, respectively, yields 

= - cos(20ci)e"^'' + KS (14) 
= vsl{'KK)dl(P + 25 sin(2(^ei)e-2^ , (15) 

which leads finally to 

2520C1 = sin(2^!)ci) cos(2(/)ei) , (16) 

where y = xj^ for the correlation length ^. The 
spaim/ dependence of cr is neglected [31]. Eq. (16) 
has the soliton array solutions [35,3] (sn, cn, dn: 
elliptic Jacobi functions) 

mi = M^{dn(ri/(fcinCm), fcin)/-R 

+ (-l)'cn(ri/(A:„,C„,),fc„,)}/2 (17) 

,5i = (-l)*^sn(n/(A;d^d),fcd) , (18) 

with = and fc,„ = fed- 

In the fits in Fig. 1 one observes that = 
^d does not hold. In addition, there is no good 
reason to neglect the spatial dependence of u 
(cf. Rcf. [36]). Considering it leads to an extra 
factor e"^'^'^ on the right side of (16) where Act 
is the difference to the fluctuations in the ground 
state. Previous results [4], (not yet self-consistent: 
fluctuations on top of the Nakano/Fukuyama so- 
lution), showed that the spatial dependence of Aa 
is indeed important. The conjecture that the nu- 
merical finding ^ ^d [3] is related to the spatial 
dependence of Act could be corroborated. 

Figs. 2 and 3 show the results of a fully self- 
consistent calculation. Clearly the spatial depen- 
dence of the fluctuations makes the kink smoother. 
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Fig. 2. Distortion of a single soliton with (solid line), with- 
out (dotted line) spatial dependence of A(7{x); fit with 
relativ width = 1.80 
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y=x/^ 

Fig. 3. Alternating magnetisation of a single soliton with 
(solid line), without (dotted line) spatial dependence of 
Acr{x); fit with relativ width = 1.45 

This was expected since the fluctuations arc max- 
imum in the center of the soliton (see Fig. 4 in 
Ref. [4]). Most interesting is the fact that the dis- 
tortion is smoothcncd more than the alternating 
magnetisation. From the fits one finds the ratio 
Cd/Cm = 1-24 in very good agreement with the nu- 
merical result at critical frustration ac (see Fig. 5 
in Ref. [3]). Hence the discrepancy between the 
existing continuum theory of spin-Peierls solitons 
and the numerical results is removed when the spa- 
tial dependence of the fluctuations are considered 
properly. 

In spite of the progress just shown the agree- 
ment to experiment is not yet perfect. The soli- 
ton widths observed experimentally still exceed the 



theoretical ones [3] . As a next step towards a quan- 
titative description we include magnetic interchain 
couplings as they are present and non-negligible in 
CuGeOs [37]. They are treated in mean-field ap- 
proximation so that only the local magnetisations 
matter. Apart from this extension the Hamiltonian 
(6) is minimised since the antiadiabatic limit leads 
to the same equations. In Fig. 4 the dependence 




0.00 0.05 0.10 0.15 0.20 



Fig. 4. Distortive (^a) a^id magnetic (^m) soliton widths at 
two average magnetisations m at a = 0.35, K = 10.2. 

of the soliton widths on an antiferromagnetic cou- 
pling Jinter between adjacent chains in a plane is 
shown. Clearly, all widths are enhanced. The ad- 
ditional coupling strongly favours the alternating 
magnetisation close to the center of the soliton. 
Thus the 2D coupling induces an increase of this 
region thereby enlarging the soliton [35]. So the 
agreement with experiment will be improved by 
taking the 2D character of the magnetic couplings 
into account. Yet further questions still remain as 
the ratio between distortive and magnetic soliton 
width is decreased on increasing Jinter in contrast 
to experiment. 

3. Summary 

In this work we contrasted the adiabatic ap- 
proach to spin-Peierls transitions with the antia- 
diabatic approach. If the latter is complemented 
by a chain-mean-field treatment both approaches 
yield very similar equations. But the interpretation 
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of the parameters differ and so do certain physical 
phenomena such as T-dependent couphngs and lo- 
cal magnetisations. For CuGeOa the antiadiabatic 
approach was advocated. 

In the incommensurate phase the form and the 
amplitudes of solitons were discussed. The exten- 
sion of the existing continuum theories as required 
by the numeric results was presented. Further steps 
towards a quantitative description of the I phase 
of CuGeOa were shown. 
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